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Running Coupling BFKL Equation and Deep Inelastic Scattering. 
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I examine the form of the solution of the BFKL equation with running couphng relevant for deep inelastic 
scattering. The evolution of structure functions is precisely determined and well described by an effective coupling 
~ l/{f3o(ln{Q^/A^) + 3.6(Qs(Q^) ln(l/a;))^'^) (until extremely small x). Corrections to the LO equation are 
relatively small, and the perturbative expansion is stable. Comparison to data via a global fit is very successful. 



Small X physics has recently been a partic- 
ular area of both experimental and theoretical 
interest. The first data on structure functions 
at very small values of x (down to x = 10~^) 
obtained by HERA [1,2] have themselves been 
enough to make this an active topic. However, the 
fact that the splitting functions and coefficient 
functions required for the construction of struc- 
ture functions have expansions containing terms 
a"(/z^)ln"'~"^(l/a;) has added extra impetus, im- 
plying that at small x one may have to account 
for leading terms in ln(l/a;) at high orders in as, 
rather than just expand naively in powers of as- 

1. The Running Coupling BFKL Equation. 

The BFKL equation [3] makes it possible to 
take account of the most leading ln(l/a::) terms 
at each order in a^, since it is the four-point 
gluon Green's function containing all the infor- 
mation about the most important small- a; be- 
haviour. At leading order, with fixed as, the 
BFKL equation is scale invariant, and has sim- 
ple eigenfunction ~ (fc^)'*', with eigenvalues x'^(7), 
where x'^il) is the Mellin transformation of the 
kernel. Transforming back to x and fc^ space the 
leading small x behaviour is driven by the saddle- 
point at 7 = 1/2, leading to small- a; behaviour of 
the form x~'^, where A = q;sx"(1/2) — as41n(2) 
{cis = Sas/ir), i.e. if as ~ 0.2, A ^ 0.5. How- 
ever, this result is potentially subject to large 
corrections. In particular the equation predicts 
large diffusion - within the gluon Green's func- 
tion the average virtuality ^ k'^ but the spread of 



important values has a width ex {as\n{l/x)y/^. 
Hence, there is considerable influence from both 
the infrared and ultraviolet regions of virtuality 
at small x. If the coupling is allowed to run this 
will clearly be very important. 

Recently the NLO BFKL equation became 
available [4,5]. Not only does this provide in- 
formation about the scale-invariant NLO correc- 
tions, but the running of the QCD coupling now 
becomes impossible to ignore - the NLO part of 
the kernel contains a term ex l3oa1[y?) ln(A:^//u^), 
where jjl is the renormalization scale, and this 
is absorbed into the renormalization group im- 
proved couphng, as{n^) — /3oQ;^(/Lt^) ln(fc^//i^) — s- 
as{k'^). Having done this, a simple way to con- 
sider the effect of the NLO corrections to the 
BFKL kernel is to consider its action on the 
LO eigenf unctions {k'^)'* ■ This results in a NLO 
eigenvalue of the form as{k'^)x'^ il) ~ ^'iif^'^)x^ il) j 
which for 7 = 1/2 leads to an intercept of 
2.8as(fc^)(l - 6.5asiP)). This is clearly a disas- 
ter, implying a very unconvergent expansion for 
the intercept. 

However, this simple approach requires modi- 
fication. The eigenvalue is not a real eigenvalue 
since it depends on k^ - the running coupling has 
broken scale invariance, changing the whole form 
of the BFKL equation. Even for the LO BFKL 
equation, the ln(fc^/A^) associated with a running 
coupling turns the 7-space BFKL equation into a 
differential equation. This can formally be solved, 
and then inverted back to fc^-space. At leading 
twist, the solution factorizes into a fc^-dependent 
part g{k'^,N) and a QQ-dependent part [6]. The 



latter is ambiguous since the 7-dependent inte- 
grand has a cut along the contour of integration 
of the inverse transformation. This is because 
the growth of the coupling in the infrared cou- 
pled with the strong infrared diffusion leads to 
an infrared divergence in the BFKL equation, and 
hence a renormalon contribution. Therefore the 
running coupling equation has no prediction for 
the input for the gluon, this being sensitive to 
nonperturbative physics. Conversely, g{k'^,N) is 
completely well-defined, being insensitive to in- 
frared diffusion. However, it is influenced the 
ultraviolet diffusion, and hence by scales greater 
than fc^, where the coupling is weak. 

Calculating dg{k'^ , N) / d\n{k'^) one extracts an 
unambiguous anomalous dimension T{N,as{k'^)) 
governing the evolution of the gluon. Consider- 
ing for the moment the use of just the LO ker- 
nel, r(iV, as(fc^)) turns out to be the usual LO 
BFKL anomalous dimension 7" (a^/A^), with a^ 
evaluated at fc^, plus an infinite series of correc- 
tions which are a power series in /3oQ:s(fc^) com- 
pared to the LO. Transforming to x space one 
finds that each term in the series behaves like 



-A(fe^) 



as X — » 0, but with accompanying pow- 



ers of as(fc^) ln(l/x) which grow as the power of 
f3oas{k'^) grows. Hence a resummation of run- 
ning coupling dependent terms is necessary. The 
series is too complicated to sum exactly, so some 
prescription must be used. I choose the ELM pre- 
scription [7], which adjusts the scale of the cou- 
pling in the LO expression to fc^ so that using 
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a,{k') = asik')-PoHkVk')aiik') + - 



in the LO splitting function one generates pre- 
cisely the NLO /3o-depcndent correction to the 
splitting function. Choosing fc^ in this man- 
ner one can find an exact expression for the 
scale choice in the effective coupling, and in the 
smafl X limit it is \n{P/A'^) = ln(k'^/A'^) + 
3.6(as(fc^)ln(l/x))i/2, (see [8] for details), which 
is entirely consistent with influence from the dif- 
fusion into the ultraviolet. The effective coupling 
is shown in flg. 1. 

It is possible to check whether this prescrip- 
tion is sensible. If the summation of the com- 
plete series is generated by this effective scale 
then the whole series is given by using as(fc^) = 
a^(fc2) - /3oln(fcVfc2)a2(fc2) _^' . . . ^^ ^^^ lO 



splitting function. This can be checked explic- 
itly for the O{{l3oas{k^))'^) term. The agree- 
ment between the term generated in this man- 
ner and that calculated explicitly is excellent un- 
til very small x indeed, i.e., the 0((/3oas(^^))^) 
term leads to a modiflcation to the scale of 
ln(fcVA2) = ln(feVA2) + 3.6(a,(fc2) ln{l/x))^/^ - 
1.2/3oas(fc2)(as(fc2)ln(l/x))2. Checks at higher 
order in /3aas{k'^) give similar results. 



Figure 1. a^-^-^ as a function of x compared to the 
constant values for two choices of ln{Q'^ /A"^). 

One can also solve for T{N,as{k'^)) numeri- 
cally and compare with the transformed split- 
ting function using the effective scale. Due to 
a zero in g{k'^,N) (flrst noticed in the context 
of a resummed kernel in [9]), the anomalous di- 
mension has a leading pole for A^ > 0, which for 
as(fc2) w 0.25 is at TV « as{k^), leading to a 
splitting function ^ gr-"s{k )^ However, the nu- 
merical T(N,as{k'^)) is in very good agreement 
with the transformation oi p^„{x,a'^'f) until one 
gets extremely near the pole - the replacement 
ln(P/A2) =ln(fcVA2) + 3.6(a,(P)ln(l/a;))i/2 is 

extremely effective until very small x indeed (see 
[10] for details). At a; ~ 10~^ the corrections 
from O{{f3oas{k'^))'^) and beyond act to stop the 
decrease of the coupling at even smaller x and 
freeze it at a^f^^^ik^) ~ OAusik^). 

The validity of using the effective coupling 



can also be checked by solving the NLO BFKL 
equation with running coupling, first defined in 
[11]. Now, as well as /3o-dependent corrections 
to the naive LO result there is also an 0{as) 
correction not involved with the running of the 
coupling, and ignoring /3o-dependent effects the 
NLO correction to the intercept of the split- 
ting function is indeed —6.5as{k^) the LO. How- 
ever, one may use the same type of prescrip- 
tion for fixing the scale in the NLO contribu- 
tion to the splitting function as at LO, finding 
that, although it is not necessary, the a'^^^ ap- 
propriate for this NLO correction is the same 
as at LO [8]. This is both for the In(fcVA^) = 
ln(fcVA2) -I- 3.6(a3(A;^)ln(l/x))i/2 behaviour for 
X > 10^^ and the freezing below this - the a; — > 
NLO corrected splitting function does not be- 
have like the naive a;-2-8"»('=')(i-6-^"»('=')), but 
like x~s.s(k')(i-2.3<iAk^)) foj. Q,^(fc2) ^ Q 25 [10]. 

Hence, the perturbative expansion is more stable. 
Exact calculations at finite x show that the NLO 
correction to the splitting function leads to fairly 
small corrections to the evolution [8]. I note that 
other authors have considered partial resumma- 
tions of the BFKL kernel which provide improved 
stability besides that associated with the running 
couphng [12,11,9]. These are necessary for single- 
scale processes, but I feel they are less important 
for structure functions than running coupling ef- 
fects. 

2. Phenomenology. 

In order to apply these results to a study of 
structure functions it is necessary to calculate the 
small X expansions for the physical splitting func- 
tions [13] which give the evolution of F2{x,Q^) 
and Fl{x, Q^) in terms of each other, in order to 
avoid factorization scheme ambiguities. This re- 
sults in al^^ of exactly the same form as above. 
In order to compare with data one can then com- 
bine these small x expansions with the normal 
expansions in powers of as in the manner in [14] . 
At present there is only sufficient information to 
work at LO in this combined expansion. How- 
ever, it is possible to use effective scales at large 
x using similar considerations as above, finding a 
coupling which grows as a; ^ 1. 



A global analysis of structure function data can 
be performed (with constraints applied for other 
data, e.g. prompt photon, high Et jets), and the 
results are very good. The x^ for 1330 data points 
is 1339 compared to 1511 for the standard NLO 
in as MRST fit [15]. Not only is the fit improved 
but, whereas the input gluon (and therefore in- 
put Fl{x,Q^)) in the conventional approach is 
valencelike at Q^ ^ IGeV^, Fl{x, Q^) in this ap- 
proach is the same general shape as F2{x, Q^) at 
this low Q^ . The prediction for Fl{x, Q^) in the 
HERA range is smaller than that using the con- 
ventional approach, being similar to [14]. 

Hence, I believe that the correct way to take 
account of the ln(l/a::) terms in the calculation of 
structure functions is to use the combined expan- 
sion for physical splitting functions, and the a'^'-' 
given by resumming /3o-dependent terms. This is 
preferred by current F2{x,Q'^) data, and predic- 
tions for other quantities are different from the 
conventional approach. 
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